Single-particle excitations and the order parameter for a trapped superfiuid Fermi gas. 
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We reveal a strong influence of a superfiuid phase transition on the character of excitations of a 
trapped neutral-atom Fermi gas. Below the transition temperature the presence of a spatially inho- 
mogeneous order parameter (gap) shifts up the excitation eigenenergies and leads to the appearance 
of in-gap excitations localized in the outer part of the gas sample. The eigenenergies become sen- 
sitive to the gas temperature and are no longer multiples of the trap frequencies. These features 
should manifest themselves in a strong change of the density oscillations induced by modulations of 
the trap frequencies and can be used for identifying the superfiuid phase transition. 



Physics of ultracold trapped atomic gases has attracted 
a lot of attention after the discovery of Bose-Einstein 
condensation As has been proven in experiments 

with Bose gases, the presence of a macroscopically oc- 
cupied quantum state (condensate) makes their behavior 
drastically different from that of ordinary gas samples. 
Trapped neutral-atom Fermi gases, being cooled to a suf- 
ficiently low temperature, should also exhibit prominent 
macroscopic quantum phenomena which are mostly re- 
lated to a superfiuid pairing phase transition: Below the 
transition point the behavior of the gas is governed by 
the presence of a macroscopic wave function of Cooper 
pairs (the order parameter). Possible versions of this 
phase transition in atomic samples have recently been 
discussed in Refs. ftfl]- However, as only a small frac- 
tion of particles is influenced by the pairing, it is not 
entirely clear how the transition will manifest itself in 
dynamic and kinetic properties of the gas. 

In this Letter we indicate a clear way of identify- 
ing the pairing transition in a trapped Fermi gas. We 
study single-particle excitations and show that they are 
strongly influenced by the pairing. Above the transition 
temperature T c the excitation eigenfrequencies are mul- 
tiples of the trap frequencies. It turns out that below T c 
the presence of the order parameter (spatially inhomo- 
geneous gap) shifts up the eigenenergies and ensures the 
existence of in-gap excitations localized in the outer part 
of the gas sample in the well formed by the gap and the 
trapping potential. Due to the temperature dependence 
of the order parameter, the eigenenergies become sensi- 
tive to the gas temperature and are no longer multiplies 
of the trap frequencies. This should strongly change the 
response of the gas to modulations of the trap frequencies 
(see below). 

The influence of pairing on single-particle excitations 
is essentially the same for all types of pairing discussed in 
Refs. 0-0], and for simplicity we confine ourselves to the 
case of the s-wave pairing. We consider a two-component 
neutral gas of fcrmionic atoms trapped in a spherically 
symmetric harmonic potential. The two (hyperfine) com- 
ponents labeled as a and are assumed to have equal 
concentrations. The Hamiltonian of the system has the 
form (h = 1) 



H = E / d ^tHo^i + V J dr^+^+if/p, (1) 

i=a,/3 

where ipi(r) with i = a, (3 are the field operators of the 
a and f3 atoms, H = — V 2 /2m + mfl 2 r 2 /2 — fi, fl is 
the trap frequency, and /i is the chemical potential which 
greatly exceeds f2 in the Thomas-Fermi limit (see e.g. 
Ref. ||) discussed below. The second term in Eq. (|l|) 
assumes an attractive elastic interaction between atoms 
in the states a and [3 (s-wave scattering length a < 0), 
V = Aira/m being the coupling constant, and m the atom 
mass. The interaction between atoms in the same hyper- 
fine states, originating in the case of fermions only from 
the scattering with orbital angular momenta Z > 1, is 
here neglected. 

The presence of a negative s-wave scattering length 
for the inter-componcnt interaction leads to a superfiuid 
phase transition via Cooper pairing in the s-wave chan- 
nel 0, with the critical temperature T c <C M- Being 
interested in the effect of this transition, we consider 
temperatures T <C A*, where the chemical potential co- 
incides with the Fermi energy in the center of the trap: 
jji ps e F = p 2 F /2m, and the Thomas- Fermi radius of 
the gas sample Rtf — Vp/Q serves as a unit of length 
{vp = pp/m). In this temperature range a small pa- 
rameter of the theory is A = 2\a\pp/ir -C I. The den- 
sity profile of the gas is n(R) — no (I — i? 2 ) 3 / 2 , where 
n o = Pp/Sir is the maximum gas density, and R the 
distance from the origin in units of Rtf- This den- 
sity profile corresponds to the local Fermi momentum 
Pf{R) = Pf{1 — R 2 ) 1 ' 2 , and the density of states on the 
local Fermi surface N(R) = mpF{R) / (2ir 2 ). 

We assume that the critical temperature T c of the 
pairing transition is much larger than Q and, hence, 
T c is very close Q to the critical temperature = 
0.28£f exp(— 1/A) in a spatially homogeneous gas with 
density no M. Below T c the gas is characterized 
by the presence of the order parameter A(R) = 
|V| (i/j a (R.)il)p(R)) . Following a standard mean-field 
procedure (see, e.g. 0), the term describing the in- 
terparticle interaction in Eq. ([!]) can be written as 
A(R)^ Q (R)^j(R) + A*(R)^(R)V£(R). Then the 
Hamiltonian (\n) becomes bilinear and can be reduced to 



1 



a diagonal form by using the Bogolyubov transformation 
generalized for the spatially inhomogeneous case fiot] : 



(MR) 
{ ^(R) 



E 



U V (B) 



I/(Ri( #t 



where a„ and f3 v are the operators of single-particle ex- 
citations. Their wave functions U V (R), V V (R.) satisfy the 
Bogolyubov-de Gennes equations 



v v 



A(R)K 
-A*{H)U V 



(2) 



with e„ > being the excitation energies. In this pa- 
per we develop a method of finding the order parameter 
A(R) and solving equations (g). 

In the vicinity of T c the order parameter can be found 
from the Ginzburg-Landau equation H : 



A{R) = 5.15 • T cV /(r c ^T)/T c exp (-R 2 /2l 



(3) 



where Z| = Ky/2\/(l + 2A) < 1, and re = 0.13(ft/T c ). 
At lower temperatures, where the Ginzburg-Landau ap- 
proach is not valid, A(R) can be obtained from the Eilen- 
berger equations [Ol which in the presence of a harmonic 
trapping potential read 



ig' u + (A • /„ - A* ■ / w )/fi = 
i/i - (2iw/ u + 2 A* ■ ffw )/n = 
if' u + (2iuf u + 2 A ■ g u )/Sl = 0, 



(4) 



with the symbol prime standing for the operator 
(pf(-R)/pf)(iVr). Eqs. (Q) follow from the well-known 
Gor'kov equations O] under the assumption ep (R) ^> T. 

The functions / w , / w are related to the anomalous Mat- 
subara Green functions in the frequency-coordinate rep- 
resentation, F = (T T V' + V' + ) an d F = (T T ipip), and the 
function g u to the combination of the normal Green func- 
tions, G = — (T T ipijj + ) — (T T ijj + ip}. In particular, 

/u>(R, n) = J d£F u (R, ( PF (R) + mt/p F (R))n) , 

where F u (TL,pii) = J dr exp (— ipnr)F w (R + r/2, R — 
r/2) is the Fourier transform of the anomalous Green 
function F w (ri,r2) with regard to r = ri — r2, and n is 
a unit vector. Similar relations are valid for f u and g u . 
The functions g u , f u and f u obey the constraint O] 



/„(R,n)/ w (R,n)-«£(R,n) = l/4 



(5) 



which also provides a normalization condition. Eqs. (^) 
should be completed by the self-consistency condition 



A(R) = \V\ -2nN(R)Tj2 J ^A.(B,n), 



(0) 



where the summation is performed over the Matsubara 
frequencies w = 7rT(2n + 1), with n being an integer. 



For real A the function f w (R., n) = / W (R, — n), and the 
solution of Eqs. (fl|) for f u , omitting the terms of order 

(ft/VT 2 + A 2 ) , can be written as 



2Vu 

+ft 2 



A 2 " + 4 

1 2(w 2 4 



ftA' 



(u 2 + A 2 ) 3 /2 

w~ ^ + A 2 )A"-5A(A') 2 
16 (w 2 + A 2 ) 7 /2 



(7) 



where the symbol prime has the same meaning as before. 
Eqs. (0) and (H) provide us with an equation for the order 
parameter A(R). Actually, the first term of Eq. (Q) 
gives a formally divergent quantity ^ w 7rT/ \Juj 2 + A 2 . 
To eliminate this divergency we renormalize the coupling 
constant V in the same way as it has been done in the 
Bogolyubov method (see, e.g., Il3)). Then we obtain 



A 



XN(R) 



A • Si/2 + 5*5/2 - 



1-i? 2 



12 



1 dA 2 - 3R 
RdR ' 1 



R 2 



'7/2- 



5 (1 -R 2 ) 
24ft 2 



fA 

dR 2 



dA 
dR 



(8) 



where S a = ttT £ w lu 2 /(lu 2 + A 2 ) a for a = 5/2,7/2 , and 

~ = i_ a r 2& 

1/2 A 7 n ^ri 0) (l- J R 2 ) io exp(f cosher)) + l' 

with 7 = 0.577 being the Euler constant. 

At T close to T c equation (8) reduces to the Ginzburg- 
Landau equation from Ref. [7]. For lower temperatures 
we solved Eq. (^J) numerically, with the boundary con- 
ditions A'(0) = and A(l) = 0. In Fig. 1 we present 
A(R) at various temperatures for A = 0.3 and = 5ft 
(T c = 0.86Tc ). The comparison with results of the 
local-density approximation shows that the latter is not 
valid in general, but gives reasonable results only at very 
low temperatures: With decreasing T, the spatial region 
where A is essentially nonzero increases and, hence, spa- 
tial derivatives of A, neglected in the local-density ap- 
proximation, become less important. 
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FIG. 1. Order parameter A(7?) at various temperatures 
r = T/T c (solid lines). The dashed lines correspond to A(R) 
in the local density approximation at r = (upper curve) and 
r = 0.99 (lower curve). 

In a spherically symmetric trapping potential elemen- 
tary excitations are characterized by the radial quantum 
number n, orbital angular momentum I, and its pro- 
jection m. The excitation wave functions can be writ- 
ten as (U V ,V V ) = i? _1 F/ m (R)(u„;(i?),w„i(i?)), and the 
Bogolyubov-de Gennes equations (0) take the form 



d 2 
dR 2 



l 



JPR 2 



Unl 
t'nl 



U n l 
-Vnl 



(9) 



where fx = (2/i/fi) > 1, £ni = £ni/Q > 0, A = A/ft, 
and the functions (u, v) are normalized by the condition 
Jq 00 {u n iu* n 'i + VniVrfi) dR = <W- At temperatures T < 
Ef ~ /i elementary excitations formed by particles with 
energies in a narrow vicinity of the Fermi surface are 
most important. In the classically accessible region of 
space the excitation wave functions exhibit strong spatial 
oscillations, with a period of order p~p^(R) <C Rtf and a 
slowly varying amplitude u n i(R), v n i(R): 



u n i 

Vnl 



exp (in pFidR 

VpMR) 



Unl 
V n l 



H.c. 



(10) 



The partial Fermi momentum is defined as pfi (R) = (1 — 
R 2 — (I + l/2) 2 /Ji 2 R 2 ) 1 / 2 , and the classically accessible 
region is specified by the condition Ri < R < R 2 , with 
the turning points i?i_2 from the equation pfi{R\,2) = 
0. Omitting the terms of order /T 1 in Eq. (|), for the 
amplitudes f± — u ± iv we obtain a pair of decoupled 
equations: 



■•I -X2 , dA 



dR 



dR 



f± = 0. (11) 



In classically inaccessible regions < R < i?i(due to 
the centrifugal barrier) and R > i?2(due to the trapping 
potential), Eqs. ( [Tl|) should be modified by replacing 
Pfi(R) with =F« \pfi{R)\, respectively, to obtain decaying 
solutions. Above T c the order parameter is zero, and 
using a standard semiclassical procedure one can obtain 
from Eqs (||)-(jll]) the well-known result = (2n + 

I + 3/2)f7 for the eigenenergies counted off the bottom 
of the potential well. In a spherical harmonic trap the 
chemical potential /i = ( j + 3/2)f2, where j is a positive 
integer. Accordingly, we obtain e^} = \2n + I — j \ Q for 
the eigenenergies of particle-like (2n + l > j, v n i = 0) and 
hole-like (2n + I < j, u n i — 0) single-particle excitations. 

Below the transition point the appearence of the or- 
der parameter A(i?) modifies the excitation spectrum. 



Just below T c the order parameter is small and exists 
only in a small spatial region of radius l& "C 1 (see Eq. 
(|)). Therefore, the presence of A(i?) only influences 
the excitations with small I and slightly shifts up their 
eigenenergies: The shifts will be of order S = A(Ri)Ia, 
i.e. much smaller than the maximum value A(0) of the 
spatially inhomogeneous gap A(R). Hence, the lowest 
excitations, namely the ones with = at T > T C1 be- 
come in-gap, i.e., have energies below the top of the gap: 
£ni ~ 5 -C A(0). With decreasing temperature, A(R) 
rapidly grows, which increases the number of in-gap ex- 
cited states. 

At temperatures well below T c the characteristic ra- 
dius of A(R) becomes of order the size of the gas sample 
Rtf, and all relevant excitations (I < m/2) are influ- 
enced by the presence of A(R). The wave functions of 
above-gap excitations (e n i > A(i?i)) extend over the en- 
tire classically accessible region R\ < R < R?. On the 
contrary, the in-gap excitations with energies e n i well be- 
low A(i?i) are essentially "expelled" from the center of 
the trap: Their wave functions are mostly localized in 
the well formed in the outer part of the sample by A(R) 
and the trapping potential (see Fig. 2). 
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FIG. 2. Wave functions u n o (solid lines) and v n o (dashed 
lines) for above-gap and in-gap excitations, obtained by nu- 
merical solution of Eqs. (8) for I = at T = 0. For illustrative 
purpose only (to reduce the number of rapid oscillations) we 
take /i = 63 instead of the actual value u ~ 970. Arrows 
indicate eigenenergies of the excitations. 

For T well below T c we have A(0) - T c /Cl > 1, and 
Eqs. ( |lT| ) for the amplitudes u n i, v n i can be solved in the 
WKB approach, omitting the terms with spatial deriva- 
tives of A(R) as they are small compared to A 2 . Then, 
in the spatial region where e > A(R) we obtain 



Unl 
Vnl 



x exp < ±i 



1/2 

T(s±u(R))- 1/2 



(12) 



PFl(R) 



dR 



Here ui(R) 



A 2 (R), and C± are numerical co- 
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efficients. In the case of above-gap excitations (e n i > 
A(iii)) these coefficients can be found by making an an- 
alytical continuation of the solution ( |l2| ) to the classi- 
cally inaccessible regions R < R\ and R > R 2 - Since for 
e n i > there are only two (classical) turning points 

i?i and i?2, the WKB quantization condition reads 



Ri 



V^-A 2 (i?) 

PFl(R) 



dR 



.(o) 



(13) 



This condition provides us with the energy spectrum of 
the above-gap excitations (see Fig. 2). The amplitudes 
u n i, v n i of their wave functions oscillate in the entire 
region Ri < R < R 2 , and decay in classically inaccessible 
regions R < Ri and R> R 2 . As well as excitations above 
T c , these excited states are doubly degenerate, but their 
energies arc shifted up by the presence of A(R). More 
accurate analysis of Eqs. (|j|) shows that the degeneracy 
is lifted, but the splitting is small. 

For the in-gap excitations the situation is more subtle 
due to the appearance of a peculiar turning point R c 
determined by the condition e n i = A(R C ). At this point 
a particle undergoes the Andreev reflection |l4|] from the 
spatially inhomogeneous gap A(i?) and transforms into 
a hole (and vice versa). This ensures a strong coupling 
between the particle-like and hole-like excitations of the 
same energy. As a result, these excitations acquire a 
superpositional particle-hole character and become non- 
degenerate, with a splitting increasing with A(0). 

In the spatial region R c < R < R 2 the amplitudes u n i , 
v n i are determined by Eq. (|l^), with the coefficients C± 
following from an analytical continuation of Eq. ( |l2| ) to 
the region R% < R < R c , In the latter region the ampli- 
tudes are given by the same Eq. ( |l2| ) where w(R) is re- 
placed by —i \uj(R)\, and the coefficients C± are obtained 
by making an analytical continuation to the classically in- 
accessible region R < R\. As a result, the quantization 
condition for the in-gap excitations, which provides us 
with their energy spectrum, reads 



(-l)^cos(20) = 2Z7(Z 4 + l), 



(14) 



where Z = \/2exp { ^A 2 (R) - e 2 nl /p Fl {R)dR}, and 

<t> = y/e 2 nl - A 2 (R)/p F i(R)dR. The wave functions 
of these excitations are mainly localized in the region 
R c < R < R 2 , where the amplitudes u„i and^ oscillate. 
For i?i < R < R c they decay exponentially. In fact, 
these amplitudes behave themselves as wave functions of 
bound states in the potential well formed by the trapping 
potential from one side and by the order parameter A(i?) 
from the other side (see Fig. 2). The existence of such 
in-gap excitations in a trapped Fermi gas originates from 
the spatial inhomogeneity of A(R). To a certain extent 
these excitations are analogous to localized states in the 
vortex core in ordinary superconductors 15|. 

For the lowest in-gap excitations the WKB approach 
for finding the eigenenergies e n i and amplitudes u n i , v n i is 



not adequate, and one has to solve Eq. ( |11[ ) numerically. 
The energies of these excitations are very sensitive to 
A(i?), and, hence, to the gas temperature. For A(R) 
in Fig. 1 at T = we find £o = 0.850 for the lowest 
excitation with I = (Eq. @ gives e = 1.06ft). With 
increasing T the value of e decreases (eo = 0.23ft for 
T = 0.99T C ) and tends to zero for T -> T c . 

In conclusion, we have found the wave functions 
and energy spectrum of single-particle excitations of a 
trapped Fermi gas and showed that they are strongly in- 
fluenced by a pairing transition. The presence of a spa- 
tially inhomogeneous gap A(R) shifts up the excitation 
energies and leads to the appearance of the in-gap exci- 
tations with wave functions expelled from the center of 
the trap. Due to the growth of the gap with decreas- 
ing temperature, the excitation eigenfrequencies become 
quite sensitive to the value of T and are no longer mul- 
tiples of the trap frequencies. These features should lead 
to a strong change of the density oscillations induced by 
modulations of the trap frequency, which can be used for 
identifying the pairing transition. For example, spheri- 
cally symmetric modulations (<5ft ~ cos(ft)) cause single- 
particle transitions between the states with the same or- 
bital angular momentum I, and above T c the amplitude 
of the density oscillations will exhibit resonances at fre- 
quencies v which are multiples of 2ft. (In this Letter we 
do not consider collective excitations related to density 
or (and) order parameter fluctuations. They are cither 
overdamped or have different frequencies.) Below T c the 
presence of A(R) changes the eigenfrequencies of some 
of the excitations, and those will not contribute to the 
density oscillations at the resonance frequencies. Hence, 
the resonance peaks broaden and become smaller. For 
A(R) in Fig. 1 already at T = 0.5T C the eigenfrequen- 
cies of all excitations with I < Jl/2 are altered by the 
pairing. Thus, the resonances in the independence of 
the density oscillations, characteristic for the gas above 
T c , will be smeared out. This feature will not be influ- 
enced by collective excitations related to fluctuations of 
the order parameter and not discussed above, since the 
eigenfrequencies of these excitations should be different 
from multiplies of 2ft. 
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